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It is shown how the initial azimuthally asymmetric flow develops from 
the free-streaming + sudden equilibration approximation to the early dy- 
namics in relativistic heavy- ion collisions. Consequences for the hydrody- 
namics and physical results are discussed. 

PACS numbers: 25.75.-q, 25.75.Dw, 25.75.Ld 

We present a description of the early-stage dynamics in relativistic heavy 
ion collisions, where the free- streaming (FS) of partons is followed by a sud- 
den equilibration (SE) to a thermalized phase, providing initial conditions 
for the subsequent hydrodynamic evolution. This FS+SE approximation 
has been proposed by Kolb, Sollfrank, and Heinz [lj. It has been further dis- 
cussed in the context of the isotropization problem by Jas and Mrowczyriski 
[2J, as well as used to analyze the early development of collective flow by 
Sinyukov, Gyulassy, Karpenko, and Nazarenko [21 El [5] • In this talk we 
point out the emergence of the initial azimuthally asymmetric flow from 
FS+SE. 

The cartoon of the approach is given in Fig. [TJ Rather than assuming 
a gradual transition from an inequilibrated partonic stage to a thermalized 
system (top panel), in FS+SE one approximates this early stage of evolution 
with free streaming followed by a sudden equilibration (bottom panel). Note 
that this is analogous to the standard treatment of the freeze-out at the 
end of the hydrodynamic evolution, where the continuous decoupling of 
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Fig. 1. Evolution of the system formed in relativistic heavy-ion collisions, consist- 
ing of partonic free streaming, equilibration, hydrodynamics, freeze-out, and free 
streaming of hadrons to detectors. Top: equilibration and freeze-out occur grad- 
ually. Bottom: the approximation of the sudden equilibration and instantaneous 
freeze-out. 



hadrons (see the talks by Yu. Sinyukov and J. Knoll in these proceedings) 
is approximated with a sudden Cooper- Frye freeze-out. 

Ever since the FS+SE has been proposed, it has been generally thought 
that it unavoidably reduces the elliptic flow, which develops hydrodynami- 
cally due to the azimuthal asymmetry of the initial density profile for non- 
central collisions. Admittedly, FS decreases the spatial azimuthal asymme- 
try of the system with time. We carefully reexamine this argument. The 
point is that while FS alone obviously cannot generate azimuthal asymmetry 
in the momentum distribution, due to the well-known fact that interactions 
among produced particles are needed to generate i?2, the SE in fact does 
the job. This is because SE is dynamical in nature, resulting in an abrupt 
change, due to interactions, of the energy-momentum tensor of the system 
into a diagonal form (in the reference frame of the fluid element) of perfect 
hydro. 

One may interpret the FS+SE approach as an approximation to vis- 
cous hydrodynamics. Instead of considering a complicated viscous system 
far from the thermal equilibrium, where the partonic cross section has a 
finite value, one initially treats the partons as free, and later supplies them 
with a large cross section which results in an instantaneous equilibration of 
the system and transition to perfect hydrodynamics. Such an interpreta- 
tion works when viscosity decreases with time, or equivalently, the partonic 
cross section increases. Confinement effects, which switch on as the dis- 
tance between the partons increases, provide a mechanism for this behavior 
and support this interpretation, although a convincing solution of the early 
thermalization or isotropization problems is still missing despite many the- 
oretical efforts. 

Here we analyze a boost-invariant system with an initial Gaussian trans- 
verse energy profile, n(xo,yo) — exp ( — 5^ ~~ ~w)-> wnere the widths a and 



b depend on centrality and are obtained with GLISSANDO [6J. We intro- 



duce the initial, to = y^o - z h an d final, r = \Jt 2 — z 2 , proper times 
of free streaming, as well as the space-time rapidities 770 = \ log fjjqpf^ and 
t\—\ log |=| . Elementary kinematics, following from the fact that a mass- 
less parton moves along a straight line with the velocity of light and a 
four momentum p^ = (pTCOshY,pT cos 0,Pt sin ^p^shmY), relates the ini- 
tial and final coordinates of the parton: 

rsinh(?7 — Y)= Tosinh(?7o —Y), x — xo + dcos</>, y = yo + dsin0, 

rf = = TCOSh(y - 7/) - ^T 2 + T 2 sinh 2 (y-7/). (l) 

Consequently, the phase-space densities of partons at the proper times To 
and r are related, 

- 1 drjodxodyo-—-; — v U J - x (2) 



dYd 2 pTdrjdxdy J dYd 2 pTdr]odxodyo 
5(770 — Y — arcsinh[ — sinh(?7 — Y)])8(x — xq — d cos cf))8(y — yo — dsin0). 

TO 

We assume for simplicity a factorized form of the initial parton distribution, 
d 6 N( T0 ) 



dYd 2 pTdr]odxodyo 



n(x 0l y )F(Y - rj 0l p T ). (3) 



In Ref. [7] we show that to a very good accuracy F{Y — 770, pr) ~ S(Y — 77). 
Then, the energy-momentum tensor of the system at 77 = is 

p2tt 

T^ v = A dcf) n (x - (r - to) cos 4>,y — (r — to) sin cp) x 

(4) 
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where A is a constant from the pt integration. Next, at each point we pass 
to the local reference frame in which the components of the energy- 
momentum tensor vanish. The four- velocity needed for the appropriate 
boost is found from the condition T^ v (x, y)u v {x, y) — e(x, y)g^ v u v {x^ y), 
where e is the energy density in the local rest frame. In the left part of 
Fig. [2] we show the profile of e together with the energy-momentum tensor 
in the local rest frame (in units of e), displayed at a few points. We note 
that it has the structure very close to the case of the perfect transverse hy- 
drodynamics of massless particles [8J , where the transverse pressure is equal 
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Fig. 2. Left: Energy-density contours with the energy- momentum tensor in the 
local rest frame (in units of e) shown at a few points (dots). Right: the profile of 

the transverse velocity, v = \/v%+v%, in units of c. 
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Fig. 3. Left: The schematic development of the partonic elliptic flow e p from 
FS+SE. Right: the value of the generated momentum asymmetry e p plotted as 
a function of the proper time when the Landau matching is imposed. The dotted 
(dashed) line corresponds to matching to isotropic (transverse) hydro, while the 
solid line shows the results of hydrodynamics only, with no FS present. The top 
curve shows the spatial asymmetry e, which decreases with time. 



to s/2. Small departures from this form, present in our case, have the struc- 
ture of the shear tensor used to include the viscosity effects in transverse 
hydrodynamics [9]. 

The right part of Fig. [2] shows the transverse velocity profile. We note 
that it is azimuthally asymmetric (stronger in the direction of the impact 
parameter), which simply reflects the original geometry. In fact, for low 
free-streaming times and close to the origin one finds v = — ~(r — ro)Vn/n. 
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Fig. 4. The pion HBT radii R S ide , ^out , ^long, and the ratio i^out/^side for central 
collisions. The darker (lighter) lines describe the results with (without) FS+SE. 
The data from Ref. [10] . 

Thus the space - collective velocity correlations are induced. 

Following Ref. [I] , we now consider a convenient measure of the momen- 
tum anisotropy, e p = ((T xx ) — (T yy )) / \{T XX ) + (T yy )), where brackets denote 
the spatial integration. In the FS phase identically e p = 0, as no interac- 
tions have occurred. Then, at the proper time r, SE occurs. As a re- 
sult, the energy momentum tensor in the local frame is replaced, due to 
interactions, with a diagonal form of the perfect hydrodynamics: — > 
diag(^, e/3, e/3, e/3) for the isotropic, or T^ u —> diag(s, e/2, e/2, 0) for the 
transverse hydro. It is this Landau matching condition which causes e p to 
jump to a nonzero value. We show the result in Fig. [31 

Finally, we compare the physical results obtained with (isotropic) hydro- 
dynamics only (starting at an early proper time To = 0.25 fm) to the results 
obtained with FS from to to r = 1.0 fm, followed by SE and (isotropic) hy- 
dro. We use the realistic equation of state [TT] . The statistical hadronization 
is carried out with the help of THERMINAT0R [T2] . In Fig. H we notice strik- 
ingly similar results for the two considered cases, not to mention the very 
good description of the HBT data, including the infamous ratio i?out/^side- 
The similarity and agreement is similar for the pr-spectra and V2 [HE], as 
well as for other centrality classes. In Ref. [13J we have also shown that our 
model calculations reproduce very well the azHBT pion interferometry. 



A practical conclusion from our study is that the inclusion of FS+SE 
may be used to delay the start of perfect hydrodynamics to "comfortable" 
times of about 1 fm/c. The physical results remain practically unaltered, 
since the decrease of the spatial anisotropy with time, resulting in milder 
hydrodynamic development of ^2, is intertwined with the buildup of the 
initial azimuthally asymmetric flow. Recall that in some studies, in order 
to obtain a proper description of the particle spectra and femtoscopy, hydro 
was used with initialization times of 0.1 fm/c [HI [15] (with viscous hydro). 
Recently, the phenomenological relevance of the initial flow, first examined 
in hydrodynamics in Ref. [16], has been emphasized in Refs. [151 E] ( see 
also the contribution of S. Pratt to these proceedings). 

Two of us (WB and WF) are grateful to Piotr Bozek and Stanislaw 
Mrowczyriski for useful conversations. 
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